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Description 

Method and electronic device used to svnthesise the sound of church 
organ flue pipes, bv taking advantage of the physical modelling 
technique of acoustic instruments. 

The present patent application refers to a method and electronic device used 
to synthesise the sound of church organ flue pipes, by taking advantage of 
the physical modelling technique of acoustic instruments. 
Numerous numerical algorithms of physical-mathematical models have been 
5 developed based on the examination of the physical behaviour of organ flue 
pipes and the sound they produce, in order to synthesise the sound emission 
of aerophone instruments in real time. 

These models are based on the mutual symbiotic interaction between a non- 
linear active section, generally defined as "excitation", and a linear passive 
10 section, generally defined as "resonator". 

The relative numerical algorithm extemporarily produces a sequence that 
represents the sound of the instrument analysed and translated into a 
physical model. 

The sound is characterised by an initial time interval, defined as "attack 
15 transient", during which intensity increases up to a certain value. 

The intensity value is indefinitely maintained over time during the second 
phase, defined as "sustain phase", during which the waveform is 
approximately periodic. 

The analytical characteristics of this waveform, of which the most important is 
20 fundamental frequency, depend on each of the parameters that regulate the 
operation of the numerical simulation. 

Being the simulation performed in the time domain instead of the frequency 
domain because of the presence of numerous non-linear functional blocks, 
the relation between the set of parameters and each spectral characteristic of 
25 the generated sequence is extremely difficult to establish a priori. 

The characteristics can be altered by changing the set of parameters, often 
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empirically, and then evaluating the effect of such a change a posteriori. 
In particular, the fundamental frequency also depends on the quantitative 
characteristics of excitation, and not only on the frequency response of the 
resonator; being the evolution of the sequence extremely chaotic during the 
5 attack transient phase, the phase of the fundamental frequency cannot be 
pre-determined once the sustain phase has been reached. 
These two peculiarities are unacceptable in high-polyphony electronic musical 
instruments, such as church organs. 

The present invention consists in an audio-digital synthesis system based on 
10 digital signal processors, which contains a programme of physical simulation 
of the sound generation of organ flue pipes. 

The programme is divided into three fundamental, conceptually independent 
sections: the first section generates the harmonic part of the sound; the 
second section generates the aleatory part of the sound; the third section 
15 processes these components by means of a transfer function with two inputs 
and one output, thus obtaining the sequence that represents the sound of the 
organ pipe. 

Because of the independence of the section that generates the harmonic part 
of the sound, the fundamental frequency and the phase of the whole 

20 waveform generated by the programme can be determined a priori. 

The numerical parameters of the simulation programme are partially 
contained in a static memory and partially obtained by processing information 
from an electronic musical keyboard and from a set of user controls in real 
time. They determine the fundamental characteristics of the generated sound, 

25 among which the main characteristics are pitch, intensity, time envelope, 
harmonic composition and aleatory component. 

For major clarity the description of the method and device according to the 
present invention continues with reference to the enclosed drawings, which 
are intended for purposes of illustration only and not in a limiting sense, 
30 whereby: 

- Figure 1 shows a realisation of a digital electronic musical instrument used 
to synthesise sounds of musical instruments by taking advantage of the 
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physical modelling technique of the invention. 

- Figure 2 shows the three fundamental functional blocks and relative 
interconnections of an audio digital synthesis programme of the sounds of 
church organ flue pipes according to the invention. 

5 - Figure 3 shows a flow chart that explains one of the three blocks of Fig. 2, 
according to which a sequence that represents the harmonic part of the 
sounds of church organ flue pipes according to the invention is generated. 

- Figure 4 shows a stable realisation of a digital harmonic oscillator with two 
status variables according to the invention. 

10 - Figure 5 shows a procedure used to generate the time variation of the 
operational frequency of the harmonic oscillator shown in Fig. 4 according to 
the invention. 

- Figure 6 shows a flow chart used to generate the aleatory component of the 
time progression of the operational frequency of the harmonic oscillator 

15 shown in Fig. 4 according to the invention. 

- Figure 7 shows an example of time envelope used in the generation of the 
sequence that represents the harmonic part of the sounds of flue pipes 
according to the invention. 

- Figure 8 shows a flow chart of a low frequency oscillator used in the 
20 generation of the sequence that represents the harmonic part of the sounds 

of flue pipes according to the invention. 

- Figure 9 shows a time progression composed of non-rectilinear sections, 
according to which the frequency of an oscillator can be changed without 
perceiving an alteration of timbre pitch according to the invention. 

25 - Figure 10 shows an algorithm for the generation of a pseudoimpulsive 
periodic sequence according to the invention. 

- Figure 1 1 shows a set of interconnected functional blocks that explains one 
of the three blocks of Fig. 2, according to which a sequence that represents 
the aleatory part of the sounds of church organ flue pipes according to the 

30 invention is generated. 

- Figure 12 shows a status device used to limit the difference between two 
consecutive samples of a sequence according to the invention. 
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- Figures 13 and 16 show an example of wave envelope used during the 
attack transient phase of the generation of sounds of flue pipes according to 
the invention. 

- Figure 14 shows a wave envelope used to generate the aleatory component 
5 of the sounds of flue pipes according to the invention. 

- Figure 1 5 shows an architecture that explains one of the three blocks of Fig. 
2, representing a mathematical model of the resonator of the church organ 
flue pipes according to the invention. 

- Figure 17 shows the mutual interaction between two functional blocks 
10 necessary for the realisation of a generic harmonic oscillator according to the 

invention. 

- Figure 18 shows an example of a pseudoimpulsive periodic waveform 
generated by the algorithm of Fig. 10, used to generate the aleatory 
component of the sounds of flue pipes according to the invention. 

15 - Figure 19 explains the operation of the status machine of Fig. 12 according 
to the invention. 

With reference to the aforementioned figures, the electronic musical 
instrument of the invention is physically composed of a set of components, 
whose type, arrangement and interconnection are shown in Fig. 1. 

20 The embodiment is shown for mere illustrative purposes, since it neither 
represents the central innovation element of the present patent nor the only 
and necessary realisation of an electronic musical instrument used to 
synthesise the sound the organ pipes by means of algorithms of physical- 
mathematical simulation. With reference to Fig. 1, the information from a 

25 musical keyboard (1 ) and a set of user controls (2) is processed by a control 
unit (3), which regulates the operation of a DSP (6) by means of a plurality of 
numerical parameters contained in a ROM (4). The DSP (6) executes the 
synthesis programme of the sound of the organ pipe in real time, upon 
management from the control unit (3), using a RAM (5) to write and read 

30 temporary data. The product of the synthesis programme is a numerical 
sequence that is suitably converted by a DAC (7) into the analogue signal 
representing the sound of the organ pipe, which can be reproduced with an 
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amplification system and a loudspeaker (8). The synthesis programme, which 
is the central innovation element of the present patent, includes three 
sections. Each section has a fundamental function in the numerical simulation 
of the sound emission of the organ pipe, as shown in Fig. 2. 

5 The block (9) generates a main harmonic sequence (10) composed of a 
series of harmonic lines, whose amplitude and frequency conveniently 
change over time. By using this sequence and taking advantage of part of the 
composition, the block (11) generates a pseudoaleatory signal that represents 
the chaotic component of the sound. The aforementioned sequences are the 

10 two input signals of the linear resonator (12) that models the frequency 
response of the resonant part of the multiple qualities of organ flue pipes, and 
whose output (13) is the sequence that represents the sound of the organ 
pipes. 

The block diagram of Fig. 3 is a detailed view of the functional blocks of the 
15 harmonic component generator (9). The oscillator (14) generates an 
approximately sinusoidal waveform (16). The fundamental frequency of the 
waveform changes over time within a range of values comprising the 
fundamental frequency of the generated musical note. The details of the 
embodiment of the oscillator and the criterion used to change frequency over 
20 time are illustrated below. 

The waveform (17) is obtained from the sequence (16) through the non-linear 
block (15): if the sequence (16) were exactly a sinusoidal sequence 

x[n] = sin[o 0 n], 

the sequence (1 7) would be 

25 y[n] = 2 (sin[<o 0 n]) 2 - 1 = -cos[2o 0 n] = sin[2a> 0 n - n/2], 

that is to say a sinusoid with double frequency than the sequence (16). 
Each of the two sequences (16) and (17) is amplified by the relevant 
multipliers (18a) and (18b), and limited by the functional blocks (19a) and 
(19b) to values within the ±CLIP1 and ±CLIP2 intervals. The outputs of the 

30 blocks (19a) and (19b) are multiplied by two sequences produced by the 
envelope generators (20a) and (20b), respectively, as illustrated below, and 
the resulting products are summed to the node (21). The sum is a sequence 
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produced by a series of linear and non-linear instantaneous operations 
performed on the waveform (16). If the waveform were exactly the sequence 
x[n], a sequence would be obtained in the node (21) whose spectrum would 
be formed by harmonic components multiple of ©o (including a>o). 

5 As illustrated below, the sequence (23) is a low frequency waveform, 

whose purpose Is the amplitude modulation of the harmonic sequence 
through the product (22). 

The element (24) is a delay line whose impulse response is the 
sequence 5.i[n - N]. Together with the products and the sum of the block 

10 (25), this element forms a linear filter whose impulse response is 

CBYP + CDELz N . 

The block (26) is a non-linear instantaneous function described by the 
following formula: 

/(x) = (x + x 0 ) - (x + x 0 ) 3 + y 0 
15 where x 0 and y 0 are independent parameters. The purpose of the block is to 
modify the mutual proportion between the amplitudes of the harmonic 
components of the sequence processed by the block. 

The block (27) is a band-pass filter, whose peak frequency corresponds to the 
fundamental frequency of the input sequence. The parameter Q of the filter is 

20 tuned up to obtain the fundamental frequency of the input harmonic sequence 
with excellent approximation. Moreover, being the phase response of the filter 
null in correspondence of the peak frequency, the phases of the fundamental 
frequency of the input and output signals of the filter are equal. This 
characteristic enables to sum the input and output sequences of the filter, with 

25 no elision effect in the fundamental frequency: the block (28) sums the 
sequences (weighing them with the parameters GAIND and GAINF), in order 
to alter the proportion in amplitude between the fundamental harmonic 
component and the group of all other harmonic components. The output of 
the block (28) is the main harmonic sequence (10). 

30 The sinusoidal oscillator (14) consists in a special embodiment of the ordinary 
harmonic oscillator with two status variables, with necessary measures to 
improve the robustness to the variation of the operational frequency in real 
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time. 

Fig. 4 shows the cycle of operations performed at each sampling interval on 

the two conveniently initialised status variables VAR1 and VAR2. 

The parameter F determines the frequency of the sinusoid produced by the 

5 status variables oscillator that is composed of the steps (29) and (31) in the 
ordinary configuration. The disadvantage of the ordinary configuration is that 
it cannot suffer variations of the parameter F in real time without altering the 
amplitude of the sinusoids described by the same variables, in function of the 
current value of the status variables. Moreover, depending on the numerical 

10 precision of the oscillator's status variables, reductions of the oscillation 
amplitude can occur even in stationary conditions. It is sufficient to amplify the 
variable VAR2 by a factor 1 +e (with s positive, but close to zero) by means of 
the step (30) and limit the width of the variable VAR1 by means of the step 
(32) to values within the interval ±1. Using these measures, the variable 

15 VAR1 describes a unitary amplitude sinusoid with excellent approximation. 
This variable is the output (16) of the block (14) of Fig. 3. The parameter F 
depends on the frequency f according to the relation 

F(f) = 2sin(7cf/f sr ) 

where f sr is the sampling frequency. The frequency f can vary in real time 
20 within an interval [f 0 - Af, f 0 + Af] sufficient to have the frequency changes 
perceived, without a collateral amplitude alteration. 

Having defined the deviation from the central frequency f 0 as 8f, this 
parameter changes in real time according to the scheme of Fig. 5. Likewise 
the signal (23), the signal (33) is a low frequency waveform whose purpose is 

25 the frequency modulation of the generated sinusoid; with the support of the 
variable VAR1 , the block (34) generates an aleatory waveform of "sample and 
hold" type, according to the scheme of Fig. 6. Ultimately, 5f varies according 
to a constant PITCH parameter (which, assuming a value in an arbitrary 
interval [1-5, 1+5], determines the fine tuning of the sinusoid) of an oscillating 

30 sequence (33) and the aleatory sequence (34). The block (34) is described in 
Fig. 6: every time the variable VAR1 passes from a negative value to a 
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positive value, the variable RNDPTCH is updated to a new value NEWRND, 
which is an aleatory variable with a probability density function uniformly 
distributed in the interval [1-8RNDP, 1+5RNDP], being 8RNDP an 
independent parameter. 

5 The two generators (20A), (20B) produce two 5-segment envelope signals, 
whose progression is generically illustrated in Fig. 7. T1...T4 are the time 
intervals in which the signal passes from level L0 to L1, from L1 to L2, from 
L2 to L3 and from L3 to zero, respectively. The generators start producing 
the respective envelope signals upon a "note on" event. Level L2 is 

10 maintained overtime for an indefinite interval SUSTAIN, whose end coincides 
with the corresponding "note off' event. Each of the two generators uses its 
own set of these 8 parameters. 

The signals (23) and (33) are produced by a "Low Frequency Oscillator" 
shown in Fig. 8. The generation method of the triangular waveform with 
15 unitary amplitude and frequency TRFREQ illustrated in the block (35) is 
implicit. The parameters TRFREQ, TRAMPL, TROFFSET, TRCOEFF1 and 
TRCOEFF2 determine the conformation of the two signals (23) and (33), 
whose common fundamental frequency is TRFREQ. In particular, the signal 

(32) is a triangular wave of average value TROFFSET and semi-amplitude 
20 TRAMPL, while the signal (33) is formed by sections of parabolas, as shown 

in Fig. 9. The relation between the values TRCOEFF1 , TRCOEFF2 and the 
independent parameter K is biunique. The special progression of the signal 

(33) is necessary to obtain a triangular frequency modulation as exactly as 
possible (ref. Fig. 5) around the nominal frequency fo. with equal progressions 

25 of the positive and negative semi-periods, if they are expressed in semitone 
cents. 

The architecture of the generator (11) of Fig. 2 is illustrated with details in 
Figures 10, 11 and 12. With reference to Figs. 3 and 10, the signal (16) 
produced by the sinusoidal oscillator (14) is amplified by a factor RTINGAIN, 
30 limited in amplitude by the block (36) to values within the interval ±1 , and then 
processed by the high-pass filter (37). Finally the non-linear block (38) cuts 
the signal's negative values. At the output of the block (38) the signal 
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(illustrated in Fig. 13) produced by the envelope generator (39) is summed 
and the result is multiplied by the parameter RTGAIN. The result RATE is a 
sequence of values used in the non-linear block (42) defined as "RATE 
LIMITER", which is part of the structure described in Fig. 11. With reference 

5 to Fig. 1 1 , the functional block (40) generates a white aleatory sequence, with 
a uniformly distributed probability density function processed by the low-pass 
filter (41 ). The obtained sequence is the input signal of the structure formed 
by the delay lines NBDL1, NBDL2, NBDL3, NBDL4, the sums NBS1, NBS2, 
NBS3, the multipliers NCGAIN, NBFBK and the non-linear block (42). The set 

10 formed by these elements, including the topology of interconnections, is 
defined as "NOISE BOX". The signal generated by the block (42), which is the 
output of the aforementioned set, is amplified by a factor NGAIN and 
multiplied by the signal produced by the envelope generator (43), whose time 
progression is illustrated in Fig. 14. The signal NOISE is the output of the 

15 generator (1 1 ) of Fig. 2. 

Fig. 12 describes the non-linear block "RATE LIMITER" (42) formed by the 
sums RLS1, RLS2, the limiter (44) and the unit delay element (45). The value 
memorised in the delay (45) is subtracted by means of the adder RLS1 from 
the input signal "IN"; the result is then limited to values within the interval 

20 ±RATE (being RATE the sequence generated by the network illustrated in 
Fig. 10), and finally summed again to the current delay value (45) at the node 
RLS2. The result "OUT" is memorised in the delay element (45) for a 
successive cycle. Fig. 13 shows the time progression of the envelope 
generated by the block (39): upon a "note on" event, starting from the level 

25 NBL0, the level NBL1 is reached in a time NBT, indefinitely sustained over 
time, also after the corresponding "note off' event. Fig. 14 shows the time 
progression of the sequence generated by the block (43): upon a "note on" 
event, the signal starts from the value NL0, reaches the value NL1 over a 
time NT1 and the level NL2 over a time NT2 sustained until the successive 

30 "note off" event. Upon this event the signal reaches the value zero in a time 
NT3. 

With reference to Fig. 11, the non-linear block "RATE LIMITER" (42) can be 
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replaced with a linear filter, whose gain has a progression described by the 
same sequence RATE generated by the architecture of Fig. 10, so that the 
structure "NOISE BOX" of Fig. 1 1 is a linear time-variant filter. 
With reference to Fig. 2, the outputs of the generators (9) and (11) are the 

5 inputs of the resonator (12) illustrated with details in Fig. 15. The functional 
blocks of the network (12) form a cycle of operations, along which a sequence 
of samples propagates for a potentially infinite time. The two contributions of 
the two generators (9) and (1 1) are added to this sequence, instant by instant 
in the sum nodes (46) and (48) nodes, respectively, to sustain the energy of 

10 the computed sequence. The structure of Figure 15 is the translation into a 
mathematical model of the resonant part of the organ flue pipe, defined as 
"pipework". In particular, the low-pass filter (47) emulates the dissipation of 
acoustic energy, with variable intensities in function of the frequency; the 
high-pass filter (49) attenuates all the frequency components lower than the 

15 fundamental frequency; by means of the product (51), the envelope generator 
(50) produces a signal that represents the time progression of the loop gain of 
the resonant system; the filter (52) alters the sequence phase, leaving its 
module unchanged; the factor TFBK (53) depends on the type of acoustic 
termination at the top of pipework; finally, the delay line BDELAY (54) 

20 considers the time needed by an acoustic pressure wave to cover the 
pipework from the base to the top and vice versa. The time progression of the 
signal produced by the envelope generator (50) is traced in Fig. 16: likewise 
the envelope of Fig. 13, upon a "note on" event, the signal passes from a 
value FBL0 to a value FBL1 in a time FBT, and then remains constant. The 

25 output sequence (13) is the signal emitted by the mathematical model of Fig. 
2 as a whole, that is to say the time representation of the sound emission of 
the organ flue pipes. 

The description continues with the original innovative characteristics of the 
audio digital synthesis technique of the sound of flue pipes. 
30 The literature on the generation of sounds of instruments with continuous 
sound emission, among which aerophone instruments, by means of the 
physical modelling technique, proposes solutions based on a mutual 
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interaction between a non-linear active part, normally defined as excitation 
(55), and a linear passive part, defined as resonator (56), according to the 
scheme of Fig. 17. In the case of aerophone instruments, the energy 
contributed to the system is in the form or sound pressure and the signal 

5 produced is the progression of the sound pressure wave irradiated by one or 
more suitable points of the resonator. The waveform p(t) is the progression of 
the air pressure that the performer (or the bellows, in the case of a church 
organ) exercises on the instrument mouthpiece. According to this progression 
and to the progression of the pressure w(t) in a suitable point inside the 

10 resonator, an oscillating acoustic pressure e(t) injected in the resonator is 
generated. Once the sustain phase has been reached, the pressure e(t) has 
the same fundamental frequency as the pressure w(t). Being linear (except 
for very special operation modes), the resonator can be described with an 
impulse response r(t), which generates the return signal w(t) and an impulse 

15 response h(t), which generates the output signal y(t). The latter is the time 
progression of the sound emission of the instrument. Being it a numerical 
simulation performed in the time domain instead of the frequency domain, the 
fundamental frequency of the oscillation on which the system stabilises, once 
the sustain phase has been reached, is extremely difficult to predict 

20 mathematically. This depends on the fact that the frequency depends on the 
time progression of the forcing signal e(t), and not only on the frequency 
values in which the amplitude spectrum of the impulse response of the 
resonator has the relative maximum values. In fact, any type of harmonic 
oscillator (electronic, mechanical, etc.) has this characteristic. With regard to 

25 wind instruments (including organ pipes), it is sufficient, for example, to 
increase the sound pressure to obtain an increase of the fundamental 
frequency of the acoustic wave, in addition to an intensity increase, although 
the characteristics of the resonant part remain unchanged. 
Another inevitable characteristic of the oscillating systems illustrated in Fig. 

30 17 is the unpredictability of the phase of the generated signal, once the 
sustain phase has been reached. Since the waveform p(t) used to stimulate 
the system is partially chaotic, and in any case it does not contain any 
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information about the phase of the stationary wave sustained by the 
resonator, the attack transient of the signal y(t) is always and unpredictably 
different. Therefore, although the waveform has always the same periodic 
time progression in sustain conditions, it is impossible to determine the 

5 evolutions that bring the system towards this progression. In quantitative 
terms, it is impossible to determine the phase of the fundamental frequency of 
any signal processed inside the stable oscillating system of Fig. 17, taking the 
instant when the stimulus p(t) starts as time origin. Together with the difficulty 
encountered in determining the fundamental frequency a priori, this is 

10 unacceptable in the field of high-polyphony electronic musical instruments, 
such as church organs. 

The synthesis system used in the present invention derives directly from the 
synthesis in the time domain described in general and is characterised by the 
total autonomy of the excitation signal from the signal produced by the 

15 resonator. In fact, the main harmonic sequence (10) extemporarily generated 
by the block (9) of Fig. 2 is the imitation, as faithful as possible, of the signal 
e(t) of the system of Fig. 17 (assuming that the latter is a good mathematical 
model of the flue pipe of a church organ), with the substantial difference that 
the fundamental frequency and the phase of this sequence, and consequently 

20 of the sequence produced by the system (13) as a whole, are perfectly 
determined a priori. 

The preparation of the numerical parameters of any oscillating system, as 
generically illustrated in Fig. 17, requires special sensitivity and skill, apart 
from the perfect knowledge of its mathematical model. This means that the 

25 good operation of the system may be impaired, and the system may become 
unstable or even inharmonious, if only one of the parameters has a value not 
included in a proper range. Moreover, different operational modes of the 
oscillator can be obtained only by acting simultaneously and with special 
attention on a plurality of parameters, with the risk of making the time 

30 evolution of one or more signals in transit along the functional blocks of the 
system uncontrollable. This makes the search for multiple sounds produced 
by this type of synthesis slow and difficult. On the contrary, a system with no 
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feedback between resonator and excitation, such as the system shown in Fig. 
2, enables to modify the numerical parameters of the three functional blocks 
(9), (1 1 ), (1 2) in a completely independent way, without impairing the good 
operation of the system as a whole. This allows obtaining a larger variety of 
5 sounds than the one obtained by means of a feedback loop system with equal 
complexity. 

The system of Fig. 3 shows a sequence of operations performed on the signal 
produced by the sinusoidal oscillator (14). The type and order of the 
operations are only one of the possible realisations used to generate a 

10 waveform sufficiently rich in harmonic components and provided with a 
suitable time evolution. In any case, some of the functional blocks of the 
system, such as the delay (24) and the non-linear function (26), derive from 
mathematical models of wind instruments known in the literature, without the 
need of using them. The originality of the system mainly consists in the 

15 adaptation of an ordinary oscillator with status variables to non-stationary 
operational conditions, by developing the functional blocks (30) and (32) of 
Fig. 4, in order to make the oscillator robust to the variations of the parameter 
F 2 of the block (29). 

With reference to Figures 5 and 6, the originality derives from the 
20 development of the generator (34) to obtain pleasing random frequency 
variations in real time. Assuming the factor (33) as constant, that is to say 
assuming the absence of the low frequency oscillation of the sequence 8f, the 
latter assumes a new random value at every period of the sinusoidal 
sequence VAR1. The result is a statistic uniformly distributed variation of the 
25 wave period, in terms of probability density function. The variation is 
perceived as a pleasant irregularity in the sound emission. Otherwise, if 8f 
assumes a new random value at every sampling instant, the length of every 
wave period will be described by a variable formed by the sum of N aleatory 
contributions, each of them provided with uniformly distributed probability 
30 density (N is the number of samples per period). In view of the Central Limit 
Theorem, the higher is N, the more the probability density function of this 
variable approaches a Gaussian function. The frequency variation would be 
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very irregular, since high frequency deviations would be obtained much more 
rarely than small deviations from the nominal frequency. This would be very 
unpleasant, since wave periods with very different length from the nominal 
length could be generated and perceived as sudden malfunctions of the 

5 generation model. 

The generator of the aleatory component (11) of Fig. 2 is completely original, 
and the embodiment of Fig. 10, 11 and 12 derives from the analysis of 
samples of sounds emitted by a large variety of organ flue pipes, and from 
some hypotheses on their operation physics. In particular, by analysing the 

10 spectrogram of the individual wave periods of a sample and using a much 
finer time resolution than a wave period, it can be noted that a large 
percentage of sound energy concentrates in a time interval much shorter than 
the period, always situated in the same position along the wave period. Such 
sound energy covers a frequency interval considerably higher than the 

15 interval covered as an average by a plurality of periods. Therefore, the 
characteristics of the spectrogram of the stationary part of the sound of flue 
pipes are similar to the spectrogram of a train of equidistant impulses, with 
the energy of the individual period concentrated in each impulse. These 
considerations justify the architecture illustrated in Fig. 10: the sequence 

20 RATE is obtained through a series of elementary deterministic operations 
performed on the sinusoid (16). Once the sustain phase has been reached, 
the sequence RATE assumes a qualitatively impulsive progression, of which 
Fig. 18 shows one example, where T 0 is the period of the sinusoid (16). 
Regardless of the method used to obtain a pseudoimpulsive sequence, the 

25 sequence is conceptually one of the inventive foundations of the generator 
(11). 

The structure illustrated in Fig. 11 is formed by the four delay lines NBDL1, 
NBDL2, NBDL3, and NBDL4. Together with the sums NBS1, NBS2, and the 
product NCGAIN, the first three delay lines form a FIR filter. The output of this 
30 filter (that is to say the sum NBS2) is processed by the non-linear element 
(42) and then, after being multiplied by NBFBK and after passing through the 
fourth delay line and the sum NBS3, injected again in the aforementioned 
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filter. If it weren't for the element (42), the structure "NOISE BOX" would be a 
linear filter, whose spectrum would have a voluntarily inharmonious 
progression, with a plurality of resonance peaks distributed in a non- 
deterministic way, depending on the length of the delay lines and the two 

5 independent parameters NCGAIN and NBFBK. These four quantities are 
dimensioned in order to imitate the frequency response of a resonator with 
irregular geometry, such as the portion of space of the organ pipe 
immediately inside the mouth. Because of the periodic oscillation of the 
sequence RATE, the element (42) causes a continuous periodic variation 

10 over time of the "gain" (not strictly definable as such, since the "RATE 
LI MITER" is a non-linear block) of the entire "NOISE BOX". In particular, with 
reference to Fig. 18, when the sequence RATE assumes the minimum value, 
the non-linear distortions caused by the block (42) imply energy losses that 
heavily reduce the resonance effects of the "NOISE BOX". Vice versa, during 

15 the (much shorter) instants in which the sequence RATE assumes relatively 
high values, the resonant effect of the "NOISE BOX" emerges and the 
intensity of the aleatory component increases. It can be noted that during the 
attack transient, because of the envelope generator (39), whose progression 
is shown in Fig. 13, the sequence RATE assumes higher values than during 

20 the sustain phase; this increases the resonance of the "NOISE BOX" during 
the first instants of synthesis, in order to simulate the acoustic phenomena 
defined as chiff, cough, etc. produced by the flue pipes if the valve that 
regulates the passage of air from the bellows to the foot is opened rapidly. 
The non linear block (42) is formed by the two adders RLS1, RLS2, the limiter 

25 (44) and the unit delay element (45). At every sampling instant, the difference 
between the previous output value and the current input value is first limited in 
width to values within the interval ±RATE and summed again to the previous 
output value, thus obtaining the current output value. The output sequence 
"follows" the input sequence, maintaining an inclination limited according to 

30 the value RATE. For mere illustrative purposes, Fig. 19 shows a chart of an 
input sequence (continuous line) and an output sequence (dotted line). In the 
instant to the inclination of the sequence IN exceeds the value RATE/sample, 
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therefore the sequence OUT separates until it re-joins at point ti, after which 
the sequence IN remains constant. In the instant t 2 the excessive inclination 
of the sequence IN causes the immediate separation of the sequence OUT 
up to the re-conjunction point U- With respect to a linear filter, the advantage 
5 of the "RATE LI MITER" is the elimination of possible discontinuities of the 
aleatory sequence, while still maintaining a sufficient bandwidth, which are 
extremely unpleasant for the human hearing. This aspect represents the 
originality of the "RATE LIMITER". 

The non linear block (42) can be replaced with any functional block whose 

10 effect on the structure "NOISE BOX" of Fig. 11 is the quantitative resonance 
variation generated by the structure, according to a periodic progression. 
As regards the linear resonator (12), the physical considerations that involve 
the choice of the functional blocks of Fig. 15 are described herein. The 
resonant part of an organ pipe, defined as pipework, can be mathematically 

15 . described, in the most elementary way, with a "comb" filter 1/(1 - FBKz~ N ), in 
which the feedback coefficient FBK is related to the loop gain of the filter and 
the parameter N is inversely proportional to the first resonance frequency of 
the same. The more complex resonator of Fig. 15 derives from this base, 
which is very used in the field of audio digital processing. Among the 

20 elements of the resonator, the function of the delay line (54) appears evident. 
The response in module of the low-pass filter (47) is designed so as to 
consider the different energy losses suffered by the various harmonic 
components during their transit along the pipework, while the high-pass filter 
(49), whose cut-off frequency is lower than the fundamental frequency of the 

25 resonator, completely eliminates the continuous component of the stationary 
wave, to take into account the fact that the average acoustic pressure inside a 
pipework is approximately equal to the external pressure. Because of the 
envelope generator (50), during the first operation phase of the resonator, the 
loop gain of the system is moderately overabundant with respect to the value 

30 once the sustain phase has been reached, in order to obtain a faster initial 
energy accumulation in the resonator, that is to say a faster attack transient of 
the generated sound. The sign of the factor TFBK (53) is especially important: 



WO 2004/086353 



17 



PCT/IT2004/000153 



a positive sign for a pipework open at the mouth and on top, and a negative 
sign for a pipework open at the mouth and closed on top. This derives from 
the physics of the reflection of an acoustic pressure wave in correspondence 
of the pipework terminations. This physical law also justifies the use of the all- 

5 pass filter (52), the most important element of the resonator from the 
conceptual point of view. If, on one side, the mono-dimensional model of the 
pipework is sufficiently accurate to justify the use of an individual delay line to 
simulate the longitudinal propagation of an acoustic wave in the pipework, the 
approximation becomes unacceptable in the wave reflection in 

10 correspondence of a structural discontinuity characterized by non-negligible 
transversal dimensions, such as the top of the pipework. The all-pass filter 
(52) modifies the total phase delay of the closed cycle formed by the 
elements (46) ... (54) in a selective way with respect to the frequency, in 
order to make the resonance of the linear resonator (12) realistically 

15 inharmonious. The same filter is optionally used to modify the value of the first 
resonance frequency of the pipework in real time through controlled variations 
of its coefficients upon a "note off' event, to simulate the phenomenon of the 
moderate reduction of the fundamental frequency of the sound of small flue 
pipes when the air inlet valve closes. 



